 
function [] = test_step_destep()

cvs = importdata('..\..\genblep\sinc-min-phase-a500-pal.csv');
d1 = cvs(:, 1)';
d2 = step_calc(d1);

wblip = get_blip('a1200', 'unfiltered', 'winuae');
wblep = get_blep(wblip);

figure(1);
hold on;
plot(wblep', '.');
plot(d2, 'g');
title('step');
legend('python blep', 'python blip - my step calc');

error_calc(wblep, d2);

d3 = destep_calc(d2);

figure(2);
hold on;
plot(d1', '.');
plot(d3, 'g');
title('destep');
legend('python blip', 'python blip - my step calc - my destep'); 

error_calc(d1, d3);

end
